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基于 人 工 神 经 网 络 模 型 的 非 球形 颗粒 忠 力 系数 预测 
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摘要 : 本 文通 过 人 工 神经 网 络 预 测 方法 对 非 球形 颗粒 气 固 上 忠 力 系数 进行 了 预测 及 分 析 。 首 先 比较 了 BP (Backpropagation) 神 
经 网 络 模型 和 RBF (Radical Basis Function) 基 神 经 网 络 模型 对 Pettyjohn 和 Christiansen 等 人 实验 工 况 中 的 结果 进行 了 预测 。 

结果 表明 ， 采 用 RBF 方法 预测 非 球形 颗粒 气 固 电力 系数 误差 较 小 ， 计 算 效率 较 高 。 同 时 ， 应 用 RBF 基 神 经 网 络 模型 ， 对 不 同 
形状 因子 下 的 气 固 忠 力 系 数 进行 了 预测 和 分 析 。 研 究 结 果 表 明 ， 人 工 神 经 网 络 可 以 用 于 非 球形 颗粒 气 固 上 忠 力 系数 的 预测 
究 ， 本 文 研究 结果 为 复杂 形状 颗粒 气 固 上 忠 力 系数 的 预测 提供 了 一 种 有 效 的 手段 。 
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Prediction of Drag Coefficient Using Artificial Neural Network in a 
Non-spher ical Particle System 
YAN Sheng-Nan TANG Tian-Qi REN An-Xing HE Yu-Rong 
(School of Energy Science and Engineering, Harbin Institute of Technology, Harbin 150001, China) 


Abstract: In this paper, the prediction and investigation on drag coefficient of non-spherical particles is presented 
applying artificial neural network. The performance between BP (Backpropagation) model and RBF (Radical Basis 
Function) model is compared. The RBF model is employed to predict the drag coefficient of non-spherical particles 
with higher efficiency and less error. The simulation results are compared with the experimental results in the 
literature by using RBF model. It reveals that artificial neural network can be applied to the prediction on the drag 
coefficient of non-spherical particles. 
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用 。 研 究 者 们 通过 构建 电力 模型 的 方法 来 描述 气 固 
0 引 襄 之 间 的 动量 交换 ， 其 中 昌 力 系数 的 选择 和 计算 对 电 
辐 体 颗粒 在 连续 性 流体 中 的 输 运 是 一 种 普遍 现 。 。 力 计 算 的 准确 与 否 有 着 不 可 忽视 的 作用 。 

象 。 针 对 固体 颗粒 的 输 运 过 程 ， 研 究 者 们 从 实验 和 根据 目前 的 研究 结果 可 以 发 现 ， 针 对 颗粒 在 流 
数值 模拟 两 方面 进行 了 大 量 的 研究 023。 但 是 ， 研 究 体 中 自由 运动 时 电力 系数 的 研究 ， 实 验方 法 是 一 种 
多 集中 于 球形 颗粒 ， 而 实际 工业 过 程 中 不 可 避免 的 。 比较 容易 的 方法 。 研 究 者 们 主要 采用 沉降 管 和 风 洞 
会 有 非 球形 颗粒 的 存在 。 由 于 非 球形 颗粒 形状 因子 。 ”实验 等 方法 对 非 球形 颗粒 的 夏 力 系数 进行 预测 世 
的 存在 ， 其 运动 特性 与 球形 颗粒 有 着 较 大 的 差异 。 。 4。 但 是 ， 采 用 实验 测量 的 方法 也 存在 不 足 ， 例 如 

办 此， 研究 非 球形 颗粒 在 连续 性 流 体 中 的 运动 状态 。 采用 风 洞 方法 时 存在 实验 空间 较 大 ， 调 试 安装 工作 
十 分 必要 。 颗 粒 在 流体 的 输 运 过 程 中 ， 气 国之 间 动 ” ” 较 复杂 ， 人 力 物力 耗费 巨大 等 。 
量 交换 对 颗粒 运动 行为 的 控制 起 着 至 关 重 要 的 作 近 些 年 来 ， 人 工 神经 网 络 模型 取得 了 长 足 的 发 
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展 ， 由 于 其 具有 能 够 精确 


型 效率 高 ， 对 于 不 完整 的 有 
效率 等 优点 ， 因 此 ， 人 工 


估算 复杂 问题 ， 
品 输 入 数据 具 
经 网 络 模 型 正 广泛 应 用 
| 以 有 效 避 免 弥 补 实验 方 


比 唯 算 模 
更 高 的 


于 不 同学 科 的 研究 之 中 ， 梧 
法 的 不 足 。 目 前 ， 在 流 态 化 领域 中 ， 应 用 神经 网 络 
模型 进行 相关 研究 已 及 


昌 关 报道 。 研 究 者 们 采用 神 


经 网 络 模型 预测 和 分 析 流 化 床 内 氨 氧 化 物 的 生成 过 


程 ， 


泄露 现象 的 检查 和 预测 以 及 
程 ， 并 取得 了 良好 的 效果 [61。 但 采用 人 工 神 


固体 垃圾 气 化 等 过 


经 网 络 


模型 对 气 固 忠 力 系 数 进行 研究 的 工作 ， 还 鲜 有 报 


、 认 
由 。 


效果 ， 
测 误差 ; 


力 系 数 进行 预测 和 对 比分 析 ， 给 出 了 不 同 


接着 ， 


综 上 所 述 ， 采 用 神经 网 络 模型 预测 非 球 形 颗 粒 
的 气 固 忠 力 系数 的 研究 还 很 少见 。 因 此 ， 本 文采 用 
了 神经 网 络 模型 预测 非 球形 颗粒 的 气 固 
首先 ， 对 比 了 BP 和 RBF 两 种 神经 网 络 模型 的 预测 
并 讨论 了 BP 和 RBF 两 种 神经 网 络 模型 的 预 
采用 RBF 模型 对 非 球 形 颗粒 气 固 忠 


电力 系数 。 


状 因子 


情况 下 ， 上 电力 系 数 随 雷诺 数 的 变化 情况 。 
1 人 工 神 经 网 络 模型 


人 工 神 经 网 络 模型 (Artificial Neural Network， 


即 ANN)， 通 常 简 
物 神经 网 络 的 启示 下 建立 的 数 志 
仿 人 脑 结构 及 其 ] 


称 为 神经 网 络 模型 ， 


人 


型 是 


种 计算 工 


*， 日 


日 石 


的 、 由 许多 碳 


究 者 定义 的 


是 一 种 在 生 


中 处 理 模型 ， 旨 在 模 
功能 的 信息 处 理 系 统 。 神 经 网 络 模 
简单 、 高 度 内 在 连接 

经 元 组 成 。 神 经 网 络 方 


法 通过 对 外 在 输入 信号 的 动力 系统 响应 ， 能 够 学 习 
高 度 非 线性 关系 与 过 程 信息 。D3 


1.1 BP 人 工 神 经 网 络 方法 


1974 年 ，Werbo 等 人 00 首 次 提 昌 
络 学 习 的 BP 算法 ， 为 多 层 神经 网 络 


实现 提供 了 一 


Rumelhart 和 M 
差 反 向 传播 算法 进行 了 详 


经 网 络 的 学 习 问 题 ， 进 一 步 ] 
展 。BP 网 络 的 拓扑 结构 包括 输入 
能 够 在 不 知道 输入 输出 具体 数学 表达 式 的 情况 
通过 学 习 来 存储 复杂 的 映射 关系 。 在 BP 网 络 
中 ， 参 数 的 学 习 通 常 采 


了 部 


上 了 用 于 神经 网 


的 学 习 训 练 与 
切实 可 行 的 解决 途径 。1986 年 ， 
cCelland 等 四 科学 家 对 多 层 网 络 的 误 
民 的 分 析 ， 解 决 了 多 层 神 
推动 了 BP 算法 的 发 


慨 、 隐 含 层 和 输出 


误差 反 向 传播 算法 ， 数 据 


自 输 入 层 经 隐 含 层 逐 层 向 后 传播 ， 网 络 的 连接 权 值 
通过 输出 层 经 过 中 间 各 层 逐 层 向 前 ， 沿 着 误差 性 能 
函数 梯度 的 反方 向 修正 ， 原 理 采 用 最 速 下 降 法 。BP 
网 络 模型 如 图 1 所 示 。 


Xl 


图 1 BP 人工 


经 网 络 模型 示意 图 nl 


Fig. 1 Schematic diagram of BP neutral network model 


BP 算法 的 提出 在 一 定 程度 上 解决 了 多 层 网 络 参 
数 训练 难 的 问题 ， 但 是 其 自身 也 有 局 限 性 。 首 先 ， 
BP 算法 需要 较 多 参数 ， 如 网 络 层 数 、 每 一 层 的 神经 
元 个 数 及 权 值 ， 且 缺乏 有 效 的 参数 选择 方法 ， 其 
次 ，BP 算法 采用 最 速 梯度 下 降 的 优化 思想 ， 在 解决 
实际 问题 时 ， 会 经 常 陷 入 到 局 部 极 小 值 中 ， 导 致 算 
法 陷入 局 部 最 优 ， 再 次 ，BP 算法 对 于 样本 的 依赖 性 
较 高 ， 如 果 样 本 集合 代表 性 差 、 存 在 大 量 矛 盾 样 本 
和 元 余 样 本 ， 网 络 的 性 能 就 会 较 差 。 最 后 ， 对 于 一 
些 复杂 网 络 的 优化 问题 ，BP 算法 受到 学 习 速 率 的 限 
制 需要 花费 几 个 小 时 。 因 此 ， 随 着 神 将 网 络 模型 的 
发 展 ， 研 究 者 们 逐渐 提出 了 新 的 神经 网 络 模型 以 弥 
补 其 不 足 。 


1. 2 RBF 基 神 经 网 络 方法 

1988 年 ， 基 于 生物 神经 元 的 局 部 响应 特性 
Broomhead 和 Lowe 将 径 向 基 函 数 引入 到 了 神经 网 
络 的 设计 中 [00， 提 出 了 径 向 基 神 经 网 络 RBF 算 
法 。 后 来 ，Jackson 和 ParkI112 分 别 于 1989 年 和 
1991 年 对 RBF 算法 在 非 线 性 连续 函数 上 的 一 致 逼 
近 性 能 进行 了 论证 。RBF 神经 网 络 是 一 种 3 层 的 前 
向 网 络 ， 其 基本 工作 原理 是 : 利用 RBF 构成 的 非 线 
性 隐 含 层 空 间 对 低 维 的 输入 矢量 进行 投影 ， 将 数据 
变换 到 高 维 空间 中 去 ， 以 使 原来 线性 不 可 分 的 问题 
能 够 变 得 线性 可 分 ， 最 终 得 到 线性 的 输出 层 。 图 2 
为 径 向 基 神 经 网 络 的 基本 结构 示意 网 。 
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图 2RBF 人 


神经 网 络 模型 示意 图 


2 


[10] 


Fig. 2 Schematic diagram of RBF neutral network model 


RBF 网 络 有 很 快 的 学 习 收 敛 速度 ， 一 个 很 重要 


的 原因 在 于 
学 习 隐 含 层 权 值 的 过 程 ， 


结构 简单 ， 有 较 快 的 收敛 速度 ， 没 有 
省 去 了 误差 在 网 络 中 逐 层 


传递 的 耗 时 过 程 。RBF 网 络 的 研究 与 应 用 ， 标 志 着 


神经 网 络 真正 走向 实用 化 过 程 ， 


目前 已 被 广泛 应 用 


于 非 线性 函数 逼近 、 模 式 分 类 、 控 制 系统 建 模 、 时 


变数 据 分 析 和 故障 分 析 诊 断 等 工程 领域 。 
本 文 的 输入 参数 为 雷诺 数 Re 和 形状 
预期 输出 值 为 气 固 忠 力 系 数 。 采 用 已 有 文 


献 中 


内 子 2， 
的 部 


分 实验 值 作为 学 习 和 训练 数据 ， 预 测 不 同 
下 气 固 遇 力 系数 随 雷 诺 数 的 变化 情况 。 


2 结果 与 讨论 


2.1 BP 与 RBF 算法 计算 结果 对 比 


状 


ROSY 


因子 


本 节 首 先 应 用 BP 和 了 RBF 两 种 方法 对 文献 中 不 


同形 状 因子 的 点 力 系数 的 预测 结果 进行 了 


Pettyjohn 和 Christiansen 等 人 [3] 的 实验 分 别 给 出 


形状 因 
电力 系数 分 布 的 实验 结果 。 本 文 将 形状 因 


对 比 。 


子 为 


0.670、0.846 和 1.000 时 


的 实验 结果 用 于 学 习 和 


了 


子 为 0.670、0.806、0.846、0.906 和 1.000 时 


训 


练 ， 通 过 训练 ， 预 测 以 上 5 种 形状 因子 下 电力 系数 


的 分 布 ， 并 与 实验 结果 进行 对 比 。 对 比 结 
子 为 0.906 时 ， 无 论 是 采用 BP 方法 


所 示 。 当 形状 因 


果 如 医 
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还 是 RBF 方法 ， 预 测 得 到 的 忠 力 系数 与 实验 结果 均 
较为 接近 ， 但 采用 BP 方法 要 比 采 用 RBF 方法 消耗 
的 时 间 较 长 ， 如 表 1 所 示 。 当 形状 因子 为 0.806 
时 ， 如 图 3 〈b) 所 示 ， 采 用 了 RBF 方法 的 预测 结果 


与 实验 结果 符合 更 好 。 


采用 BP 方法 和 RBF 方法 对 不 同形 状 因 


力 系数 预测 的 计算 时 间 对 比如 表 1 所 示 。 


村 
从 表 | 


时 碟 
Ph 可 


以 明显 看 出 ， 在 相同 运行 环境 下 ， 采 用 BP 方法 的 


计算 时 间 平 均 在 80s 左右 ， 采 用 RBF 方法 的 计算 时 


间 在 8s 左 
相同 数据 量 作 为 学 习 和 


四 


右 ， 是 BP 方法 的 10%。 由 


表 1 BP 和 RBF 计算 时 间 对 比 


Table 1 Comparison of computing time 


此 可 见 ， 采 用 
| 练 数 据 预 测 相同 形状 因子 
情况 下 的 忠 力 系数 时 ，RBF 方法 的 计算 效率 更 高 。 


形状 因子 1.000 0.906 0.846 0.806 0.670 
BP 83s 87s 85s 84s 84s 
RBF 8.41s 8.93s 9.01s 8.90s 8.14s 
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果 与 实验 结果 的 平均 误差 分 布 情 
] RBF 方法 的 预测 误差 始终 小 于 10%， 
] BP 方法 预测 结果 


显 看 出 ， 采 
预测 误差 分 布 较为 稳定 。 而 采 


(b) $ =0.806 


图 3 BP 与 RBF 方法 预测 结 
Fig. 3 Comparison of simulation results by BP and RBF models 


图 4 给 出 了 采用 不 同形 状 因 


子 颗粒 时 ， 预 测 结 


况 


果 对 比 


。 从 图 中 可 以 


误差 分 布 波动 较 大 ， 特 别 是 当 形 


时 ， 


状 


因子 为 0.670 


明 


BP 方法 的 预测 误差 接近 21%， 预 测 误差 分 布 
介 于 10% 和 20% 之 间 ， 误 差 较 大 ， 预 测 结 果 的 可 信 
性 较 低 。 


综 上 所 述 ， 采 用 BRF 方法 预测 结果 更 接近 于 实 
验 值 ， 预 测 误差 小 ， 效 率 高 ， 对 于 不 同形 状 因子 颗 
粒 气 固 忠 力 系 数 的 预测 是 一 种 较为 合理 且 有 效 的 方 
法 ， 可 以 为 不 同形 状 因子 情况 下 的 点 力 系 数 的 选取 
提供 更 合理 的 理论 依据 。 因 此 ， 本 文 将 采用 RBF 的 
方法 对 不 同形 状 因子 颗粒 气 固 息 力 系数 进行 预测 和 
分 析 。 
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图 4 BP 与 RBF 预测 误差 分 析 

Fig. 4 Error analysis of BP and RBF models 
2. 2 不 同形 状 因子 遇 力 系数 预测 
由 上 文 对 比 可 知 ，RBF 方法 可 以 很 好 的 预测 不 
同形 状 因子 情况 下 的 气 固 上 忠 力 系数 ， 并 且 耗 时 较 
短 。Pettyjohn 和 Christiansen 等 人 的 实验 结果 共 提 
供 了 5 组 实验 数据 ， 本 文选 取 其 中 3 组 用 于 训练 和 
学 习 ， 并 预测 用 于 训练 和 学 习 的 形状 因子 情况 下 的 
电力 系数 ， 与 实验 结果 进行 了 对 比 ， 对 比 结果 如 图 
5 所 示 。 从 图 中 可 以 发 现 ， 当 形状 因子 为 1 时 预测 
结果 与 实验 结果 对 比较 好 ， 说 明 对 用 于 训练 和 学 习 
的 数据 进行 预测 可 以 得 到 理想 的 预测 结果 。 
由 图 可 见 ， 上 忠 力 系数 随 着 雷诺 数 的 增加 而 逐渐 
降低 ， 当 雷诺 数 增加 到 100 左右 时 ， 暇 力 系数 值 基 
本 保持 不 变 ， 但 是 当 雷 诺 数 达 1000 时 ， 忠 力 系 数 
种 有 增加 。 预 测 结果 与 实验 结果 平均 误差 为 8%， 
对 比 结果 较 好 。 
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图 5 RBF 网 络 预测 暇 力 系数 与 文献 实验 数据 比较 〈 学 习 数 
据 ) 


Fig. 5 Comparisons between RBF Simulation results and 
experimental data (test data ) 


为 了 进一步 说 明 RBF 预测 方法 的 合理 性 ， 本 文 
将 训练 学 习 到 的 模型 对 研究 工 况 中 的 另外 2 组 形状 
因子 下 昌 力 系数 进行 了 预测 ， 并 与 实验 结果 进行 了 
对 比 ， 对 比 结果 如 图 6 所 示 。 从 图 中 可 以 明显 发 
现 ， 预 测 结果 与 实验 结果 吻合 较 好 。 在 雷诺 数 转折 
位 置 CRe=100) 处 ， 预 测 结果 与 实验 结果 略 有 偏 


差 ， 但 整体 预测 效果 与 实验 结果 较为 接近 ， 平 均 预 


测 误差 约 为 7% 和 8%。 
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图 6 RBF 网 络 预 测 忠 力 系 数 与 文献 实验 数据 比较 (预测 数 
据 ) 


Fig. 6 Comparisons between RBF Simulation results and 


experimental data (predicted data ) 


综 上 所 述 ，RBF 预测 结果 较为 准确 ， 可 以 为 不 
同形 状 因子 情况 下 的 忠 力 系数 进行 预测 ， 并 为 非 球 
形 颗粒 忠 力 系数 的 计算 提供 了 合理 的 理论 依据 。 

图 7 给 出 了 采用 RBF 模型 对 不 同形 状 因子 颗粒 
气 固 忠 力 系 数 进行 预测 的 结果 。 文 中 分 别 预 测 了 


SS 


状 因子 为 0.5、0.6、0.7、0.8、0.9 及 1.0 时 电力 系 
数 的 分 布 情况 。 从 图 中 可 以 看 出 ， 在 不 同形 状 因子 
下 ， 忠 力 系 数 随 雷 诺 数 的 变化 情况 相似 : Re<10 

时 ， 上 忠 力 系数 与 雷诺 数 基 本 旦 线性 变化 ， 随 着 雷诺 
数 的 增加 电力 系数 逐渐 减 小 ， 当 10<Re<100 时 ， 电 
力 系 数 变化 趋势 逐渐 趋 于 平缓 ; 当 Re>100 时 ， 电 
力 系数 基本 保持 不 变 。 对 比 不 同形 状 因子 下 电力 系 
数 的 分 布 可 以 发 现 ， 当 Re<1 时 ， 电 力 系数 彼此 较 
为 接近 ， 当 Re>1 时 ， 不 同形 状 因 子 下 ， 忠 力 系 数 
的 差距 逐渐 变 得 明显 ， 并 且 随 着 雷诺 数 的 增加 形状 
因子 对 电力 系数 的 影响 逐渐 增强 。 当 Re>100 时 ， 

不 同形 状 因 子 下 忠 力 系数 随 雷 诺 数 的 变化 逐渐 减 

小 ， 数 值 上 略 有 增加 ， 但 基本 保持 不 变 。 
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图 7 RBF 神经 网 络 方法 预测 忠 力 系数 随 形 状 因子 的 变化 情况 
Fig. 7 Influence of shape factor on drag coefficient by RBF 


model 


3 结论 


昱 


本 文采 用 人 工 神经 网 络 模型 ， 对 非 球形 颗粒 气 
固 提 力 系数 进行 了 预测 及 分 析 ， 并 将 模拟 结果 同文 
献 中 实验 数据 进行 对 比分 析 。 结 果 表 明 ， 人 工 神经 
RoR 以 用 于 非 球形 颗粒 气 固 昌 力 系数 的 预测 研 
究 。 主 要 结论 如 下 : 

(1) BP 模型 和 RBF 模型 预测 结果 与 实验 结果 
对 比 发 现 ， 采 用 RBF 模型 的 预测 结果 与 实验 结果 更 
加 接近 ， 误 差 较 小 ， 计 算 效 率 较 高 ; 

(2) 采用 RBF 方法 预测 不 同形 状 因子 下 电力 
系数 的 分 布 ， 可 以 发 现 ， 不 同形 状 因子 下 电力 系数 
随 雷诺 数 的 变化 趋势 基本 相同 ; 

(3) 在 相同 雷诺 数 下 ， 随 着 形状 因子 的 减 小 
忠 力 系数 逐渐 增 大 。 


i 


a 


综 上 所 述 ， 采 用 神经 网 络 的 方法 可 以 较 好 的 预 
测 不 同形 状 因子 下 的 气 固 忠 力 系 数 ， 同 时 预测 结果 
表明 形状 因子 对 上 忠 力 系数 的 影响 不 可 被 忽略 ， 本 文 
的 工作 对 非 球 形 颗粒 气 固 忠 力 模 型 的 预测 和 构建 具 
有 一 定 的 参考 价值 。 
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